Reconciliation of stability analysis with simulation in a bus route model 



Scott A. Hill 

James Franck Institute and Department of Physics, University of Chicago, Chicago, Illinois 60637. 

(Dated: February 1, 2008) 

In a recent paper (Physica A 296, 320 (2001)), T. Nagatani proposes a time-headway model 
for buses on a bus route, and studies the stability of that model's homogeneous solutions. While 
investigating the phase diagram gotten by varying the rate of passenger arrival and the initial spacing 
of buses, he discovers a discrepancy between the results of simulation and those of linear stability 
analysis. In this paper, we reconcile this discrepancy by noting the existence of three types of phase 
diagrams, and present simulational results which confirm the stability analysis. 



I. INTRODUCTION 

While there has been much interest in the study of automobile traffic jjj among scientists, there have been few 
corresponding studies of buses |^, ||, |J. The dynamics of a bus route, while having some similarities with those of 
general traffic, differ due to the added interaction of buses with passengers at designated bus stops. A good reason 
for studying the dynamics of bus routes is that they are so often unstable. Buses are initially spaced at regular 
intervals. However, if one bus is delayed for some reason, it will then find a larger number of passengers waiting for 
it at subsequent stops, thus delaying it further. Meanwhile, the bus following finds fewer passengers waiting for it, 
allowing it to go faster until eventually it meets up with the delayed bus. Clusters of three, four, or more buses have 
been known to form due to this dynamic, resulting in slower and less frequent service. 

In reference , Nagatani presents a time-headway model for buses which is an extension of earlier work . Using 
linear stability analysis, he is able to determine the regions of parameter space in which the homogeneous solution 
(i.e., with buses spaced evenly apart) is unstable. However, his analysis only agrees qualitatively with results found 
by direct simulation of the model. 

In this paper, we make a slight modification to Nagatani's model, and discover that one can reconcile the simulation 
results with that of stability analysis. On doing the analysis, one finds that there are three different types of phase 
diagrams possible: one which agrees with Nagatani's theoretical predictions, a second which agrees with his simulation 
results, and a third which he does not see. We present our own simulational data for these three separate cases, showing 
quantitative agreement. 



II. MODEL 

We consider the following model [|| of buses on a bus route. Bus stops are labelled by s = 1, 2, . . . where stops s 
and s + 1 are a distance L apart. There are J buses, j — 1, . . . , J, which travel from stop to stop, with bus j = J in 
the lead and bus j — 1 in the rear. Every bus stops at every stop, and buses do not pass one another. The arrival 
time tj_ s of bus j at stop s is given by 

t jiS = tj.s-i + — - h A7A*,- iS _i, (1) 

where Aij, s _i = £j,a— l — ij+i.s-i is the time-headway, the time gap in front of bus j at stop s — 1. The penultimate 
term in Eq. |l| is the time it takes for the bus to travel from stop to stop; Vj. s _i is the average velocity of the bus as 
it moves. The last term is the time it takes for passengers to board the bus at stop s — 1: A is the rate of passengers 
arriving at the bus stop, so XAt 3 ; jS _i is the number of passengers that have arrived since bus j + 1 left the stop. The 
parameter 7 is the time it takes one passenger to board the bus, and so X^fAtj ■ s -i is the amount of time needed to 
board all the passengers. One could also account for the time it takes riding passengers to leave the bus; Nagatani 
shows, however, that this is not an important effect, so we will ignore it in our own analysis. 

It is reasonable to assume that a bus driver will try to prevent bunching by slowing down when the gap between 
his bus and the next is too small. One can model[^) this discretion by writing the average speed Vj, s as a function 
V(Atj :S ) of the headway, where 

T/ , A ,, tanh[Af-* e ] +tanhf c 

V (Ai) = U roitl + (U roax - UnunJ — — — — . (2) 

1 + tann t c 

The hyperbolic tangent factor acts as a spread-out step function centered at t Cl which is roughly the minimum time 
gap that a driver prefers to have in front of her. We will set t c = 2 in all that follows. The parameter v n 
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speed the bus would travel with an infinite time headway. The parameter v m i n is currently absent from Nagatani's 
formulation although he did use it in an earlier paper [EL We reintroduce it here for generality. 

It is convenient to work, not with the arrival times tj t8 , but with the time headways At^,,. We thus rewrite Eq. |l| 
in terms of the headways: 
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III. STABILITY ANALYSIS 



We consider the stability of a homogeneous flow of buses; that is, a situation where all buses have the same 
headway Ato. NagataniQ expands Eq. || in terms of small deviations from homogeneous flow, and finds that the 
stability condition for small disturbances at long wavelengths is 

A7 Z(^)! <F - (a(0)<(A7+1) a^)!, (4) 

where Ato is the initial, constant spacing between buses, V(At) is the velocity function in Eq. || and V'(Ato) is the 
derivative of that function. We can rewrite this in terms of the passenger arrival rate A7, 

T nA*o) ^ x „ T V'^to) ( ,, 



where, specifically, 



LV'(Ato) = L{ Vmax - v mi n)(l ~ tanht c )(l - tanh 2 At ) 
^(Ato) 2 [f m in(l — tanh Ato) + u ma x(l — tanht c ) tanh Ato] 2 



(6) 



This inequality will allow us to construct a phase diagram (Fig. similar to that of Fig. 8 in Ref. 0], where we vary 
the loading rate A7 and the initial spacing between buses Ato ■ 

Depending on the choice of parameters, there are actually three general forms this phase diagram can take. The 
first, corresponding to Nagatani's analytical results, occurs when the lower bound on A7 lies entirely below the Ato- 
axis, making it irrelevant since the loading rate A7 is always non-negative. After a straightforward calculation, one 
can show that the lower boundary curve reaches its maximum for that value Ato which satisfies 

Wmin(l - tanh At ) = v max (l - tanht c ). (7) 

From this wc find that the lower boundary never rises above the axis when 
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t (l - tanht c )] 2 > L(v m&x - iw)(l - tanht c ) (8) 



This type of phase diagram is seen in Figure [l]a. When Eq. ^| does not hold, then the lower bound in Eq. |^ rises above 
the Ato axis, and one gets a phase diagram more like that in Figure ^b. 

Notice from Eq. ^ that the boundary curve has no maximum when v m i n = 0. Figure ^c shows the phase diagram 
for this case. The similarity between Fig. |l|c and the simulation results in Fig. 3 of Ref. |Q] is not unexpected, since 
Nagatani does not use a minimum velocity in this paper. The puzzle is why the phase diagram of his stability analysis 
is more similar to Figure 0a; perhaps he used v min 7^ in his analysis. 



IV. SIMULATION 



To test our theory, we evaluate Eq. £5] iteratively in s. Our initial conditions are 

At jfi = At ±0.lrj (9) 

where Tj are random numbers chosen between —1 and 1. We use periodic boundary conditions in bus number; so 
for example Atjv, s = fjv> — ti.s- By stop s = 1000, our simulations each have one of four different outcomes: each 
outcome is represented by a different symbol in Fig. |l|, and an example of each can be found in Fig. ||. Many choices 
of parameters result in two or more buses being clumped together (At^ = for at least one value of i); these are 
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FIG. 1: Phase diagrams for three different sets of parameter values. The curves are given by L y^^l and L y^^l — 1, 
which theory predicts will bound the stable regions of phase space. The symbols represent simulation results: the small dots 
(•) label those simulation runs which are stable. The crosses (x) are those runs which end up with one or more headways equal 
to zero. The stars (*) mark those runs which end in a "kink-jam" phase. The boxes (□), which only appear on the left-hand 
side of Figure b, mark those runs whose headways are roughly homogeneous but higher than the initial headway Ato. (See 
Fig. ^ for examples of these four possible outcomes.) 



labelled by crosses in the figure. Other runs result in what is known as a "kink-jam" phase y], where some of the 
headways are much larger than their initial values, while others are smaller. In a few cases, all of the headways are 
roughly homogeneous, but with an average value larger than the initial value; this we also consider unstable. When 
all of the headways remain roughly where they began, we consider the run to be stable. We see in Fig. [I] very good 
agreement between theory and simulation in all three sets of parameters. 

In summary, when one takes into consideration the minimum velocity v m i n and the fact that there are three possible 
phase diagrams for this bus system, one is able to reconcile Nagatani's simulations with his stability analysis to a 
high degree. 

The author would like to thank Dr. Raghuveer Parthasarathy for indirectly encouraging this work. This work was 
supported by the Materials Research Science and Engineering Center through Grant No. NSF DMR 9808595. 
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FIG. 2: Examples of possible outcomes of our simulation. The plot (a) shows a typical stable result, with some fluctuation 
about the initial value 1.5 (in this case). Plot (b) shows an example of those runs with many buses bunched together; note that 
this particular (not uncommon) result is completely unphysical, and arises because we have gone far from equilibrium. Plot (c) 
shows a typical "kink-jam" phase; such phases are marked by having at least one gradient \Atj, B — At, iS _i| greater than half 
the distance between maximum and minimum. The last plot, plot (d), looks similar to the stable case, but the headways have 
all drifted to be much larger than the initial value of 0.5. To be specific, we classified such runs as "unstable" if no headway 
came within 0.1 of the initial headway Aio- 
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